Unambiguous evidence for extended s-wave pairing symmetry in hole-doped 

high-temperature superconductors 



o 
o 



in 



i 

CO 

c3 



I 

C 

o 
o 



> 
m 
oo 

o 
m 
o 

03 



c 

o 
o 



X 
S3 



Guo-meng Zhao* 

Department of Physics and Astronomy, California State University at Los Angeles, Los Angeles, CA 90032, USA 

We have analyzed data from angle resolved photoemission spectroscopy, Fourier transform scan- 
ning tunnelling spectroscopy, and low-temperature thermal conductivity for optimally doped 
Bi2Sr2CaCu2 08+;/ in order to discriminate between d-wave and extended s-wave pairing symmetry. 
The combined data are inconsistent with d-wave symmetry, but quantitatively consistent with ex- 
tended s-wave symmetry with eight line nodes. We also consistently explain all the phase-sensitive 
experiments. 



The phenomenon of superconductivity involves the 
pairing of electrons into Cooper pairs [1]. The internal 
wavefunction (gap function) of these Cooper pair obey a 
certain symmetry which reflects the underlying pairing 
mechanism. It is known that conventional superconduc- 
tors (e.g., Pb and Nb) possess an s-wave gap symmetry 
that reflects the phonon mediated electron-electron pair- 
ing [1]. On the other hand, the gap symmetry of high- 
temperature cuprate superconductors has been a topic 
of intense debate for over fifteen years. Three symmetry 
contenders have been isotropic s-wave, d-wave and ex- 
tended s-wave, as depicted in Fig. 1. Both d-wave and 
extended s-wave have line nodes and change sign when 
a node is crossed. A majority of experiments testing the 
symmetry (e.g., penetration depth, thermal conductiv- 
ity, and specific heat measurements) have pointed to the 
existence of line nodes in the gap function [2-5] . Qualita- 
tively, these experiments are consistent with both d-wave 
and extended s-wave gap functions. 
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FIG. 1. Three allowed pairing symmetries appropriate for 
Cu02 planes in the high-T c superconductors. Both d-wave 
and extended s-wave have line nodes and change signs when 
a node is crossed. 

There has been much experimental evidence for a 
d-wave symmetry of superconducting condensate (or- 
der parameter) for hole-doped cuprate superconduc- 
tors. In particular, phase-sensitive experiments based 
on planar Josephson tunneling [6] appear to pro- 
vide compelling evidence for a d-wave order parame- 
ter symmetry. Further, some angle resolved photoemis- 
sion spectroscopy (ARPES) studies on nearly optimally 
doped Bi 2 Sr 2 CaCu 2 8+ y (BSCCO) [7,8] show a very 



anisotropic gap that may be consistent with d-wave sym- 
metry. As a result, there is a widespread belief that the d- 
wave gap symmetry is now firmly established. Neverthe- 
less, there is also overwhelming evidence favoring an ex- 
tended s-wave gap {A\ g symmetry) [9,10]. This evidence 
includes data from phase-sensitive experiments based on 
out-of-plane Josephson tunneling [11-13], ARPES stud- 
ies on heavily overdoped BSCCO [14], single-particle tun- 
neling spectroscopy [15], Raman spectroscopy of heavily 
overdoped cuprates [16], Nonlinear Meissner effect [17], 
and inelastic neutron scattering [18]. The measurements 
of the physical properties that are related to low en- 
ergy quasiparticlc excitations [2-5] have definitively ex- 
cluded a nodeless s-wave gap symmetry, but cannot dis- 
tinguish between d-wave and extended s-wave unless one 
makes quantitative comparisons between theory and ex- 
periment. 

Although there are more experiments favoring ex- 
tended s-wave than d-wave gap symmetry [9,10], the ex- 
tended s-wave evidence has been less well known and un- 
dervalued. If the gap symmetry for hole-doped cuprate 
superconductors is extended s-wave, can we consistently 
explain all the phase-sensitive experiments? 

The next question is: If the intrinsic bulk gap sym- 
metry is extended s-wave, can we definitively explain the 
ARPES data for BSCCO? The ARPES data reported 
in the 1995 paper [19] for a slightly overdoped BSCCO 
sample with T c = 87 K are consistent with an extended s- 
wave gap symmetry. In contrast, the other ARPES data 
reported in the 1996 paper [8] for the same sample, which 
were measured by the same group with the same energy 
and momentum resolutions, are consistent with a sim- 
ple d-wave gap function. Now a question arises: Which 
ARPES data are more reliable. One can easily check 
that both the Fermi-surface (FS) and superconducting 
gaps reported in the 1995 paper are in perfect agreement 
with those determined by independent and very precise 
Fourier transform scanning tunnelling spectroscopic (FT- 
STS) data [20] (also see Fig. 2 below). The FS reported in 
the 1995 paper is also in excellent agreement with that 
recently determined using the symmetrization method, 
the most reliable method to extract the FS from ARPES 
data [21]. This indicates that the 1995 ARPES data are 
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reliable because they are in perfect agreement with three 
independent sets of data, which are very precise and in 
perfect agreement with each other. In contrast, the FS 
and superconducting gaps reported in the 1996 paper [8] 
are well off from those determined from the very pre- 
cise FT-STS data (see Ref. [20]). Moreover, the FS re- 
ported in the 1996 paper [8] disagrees significantly with 
that determined recently from the reliable symmetriza- 
tion method [21]. This implies that the 1996 ARPES 
data are not reliable because they are not consistent with 
other three independent sets of data, which are very pre- 
cise and in perfect agreement with each other. 
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FIG. 2. The angle dependence of the superconduct- 
ing gap A(9) in the Y quadrant for slightly overdoped 
Bi 2 Sr 2 CaCu2 08+j / crystals with T c = 86-87 K. The gaps 
are extracted from ARPES data [19] (solid circles) or from 
FT-STS data [20] (open circles). The solid line is the best 
fitted curve by Eq. 1. Here 6 is the angle measured from the 
Cu-O bonding direction. 

In Fig. 2 we plot the angle dependence of the supercon- 
ducting gap A(6>) in Y-quadrant for a slightly overdoped 
BSCCO with T c = 86-87 K. The gaps are independently 
determined from ARPES [19] and FT-STS studies [20] on 
similar BSCCO crystals. It is striking that two indepen- 
dent data sets are right on top of each other in the angle 
range accessible to both ARPES and FT-STS. The other 
two sets of the ARPES data [8,22] are well off from the 
FT-STS data points [20] due to the unreliable extraction 
of the FS (see the above discussion). Only the gap near 
the antinodal direction is found to be the same in all the 
three measurements [8,19,22], suggesting that the error in 
the FS has little effect on the accuracy of the gap extrac- 
tion along the antinodal direction. Based on the reliable 
1995 ARPES data (see the above discussion), it is appar- 
ent that the gap at 9 = 45° is finite rather than zero for 
this nearly optimally doped BSCCO. This is further sup- 
ported by another set of the ARPES data for a heavily 



overdoped BSCCO with T c = 60 K, which clearly shows 
that the gap at 6 = 45° is about 9±2 mcV (Ref. [14]). Re- 
cent ARPES data for Pb-doped BSCCO with T c = 70 K 
also indicate that the kink energies in the electronic dis- 
persion along the "nodal" and antinodal directions differ 
by less than 6 meV (see Fig. 4b in Ref. [23] ) . This implies 
that the gaps along the "nodal" and antinodal directions 
differ by less than 6 meV, and that the gap along the 
"nodal" direction is about 10 meV if the gap along the 
antinodal direction is about 16 meV. These ARPES data 
thus suggest that the anisotropy between the gap at 9 
= 45° and at 9 =0 decreases with doping. This conclu- 
sion is also consistent with Raman scattering data, which 
indicate that the gap difference between 9 = and 45° 
becomes negligibly small for heavily overdoped BSCCO 
and Tl 2 Ba 2 Cu0 6+!y (Ref. [16]). 

We can fit the combined data points in Fig. 2 by an 
extended s-wave gap function: 

|A(0)| = |A(cos40 + s) + Acos80|. (1) 

Here we include the next harmonic term cos 89 of the 
extended s-wave symmetry to account for the high har- 
monic correction to the Fermi surface. One can clearly 
see that the fit is excellent with the fitting parameters: 
A =19.43(46) mcV, s = 0.874(22) and A = -2.01(41) 
meV. From the fitted curve, we find that the line nodes 
are located at 9 n = 36.7° and 53.3° in the first quadrant, 
and the maximum gap is 34.3 meV. Since the intrinsic 
bulk maximum gap for the optimally doped BSCCO is 
33-34 meV, as seen from intrinsic tunneling spectroscopy 
[24], we conclude that the top Cu02 layer of this slightly 
overdoped BSCCO is optimally doped. 

Because the ARPES determined gap near the antin- 
odal direction is very accurate with an uncertainty of 
±1.5 mcV (Ref. [19]), we can fit only the FT-STS data 
with a constraint that the fitted curve is exactly through 
one ARPES data point near the antinodal direction (the 
first ARPES data point in Fig. 2). Such a fit leads to al- 
most the same fitting parameters as the fit that includes 
all the ARPES data points. 

The FT-STS data points alone could be also consistent 
with a d-wave gap function: 

|A(0)| = |A M [(1 - 5) cos 2(9 + 5 cos 6(9] | . (2) 

Here cos 69 is the next harmonic term of the d-wave gap 
function. McElroy et al. [20] fit their FT-STS data by 
Eq. 2. The best fit leads to the fitting parameters: Am = 
39.3 meV and B = 0.182. Although the d-wave fit is also 
good, the fitted maximum gap (39.3 meV) is well above 
the value (34 ±1.5 mcV) found from ARPES (Ref. [19]). 
Because the upper limit of the antinodal gap is bounded 
by the peak position of the spectrum, which is 35 meV 
below the Femi level [19], the fitted maximum gap (39.3 
meV) is unphysical. 

In order to definitively exclude the d-wave gap symme- 
try, one could quantitatively compare both d-wave and 



2 



extended s-wave predictions with other experiments. One 
such prediction is the linear slope of the low-temperature 
electronic thermal conductivity, k /T, which is directly 
related to the Fermi velocity vf and momentum kp in the 
nodal directions, and the slope S — dA(6)/dO at nodes. 
The former two quantities can be obtained from ARPES 
data while the latter one can be readily calculated from 
the gap function deduced from a fit. The residual thermal 
conduction is due to a fluid of zero-energy quasiparticles 
induced by the pair-breaking effect of impurity scatter- 
ing near the nodes in the gap. Calculations for the heat 
transport by nodal quasiparticles in two dimensions give 
a general expression [25] 



1l = n^--(— + 

T ZTid K v 2 v F h 



(3) 



where n/d is the stacking density of Cu02 planes, v 2 = 
S/(hkp) 1 and N is the number of nodes per quadrant. 
Remarkably, theory [25] predicts that k /T is indepen- 
dent of impurity concentration and that Eq. 3 is still 
valid even when vertex and Fermi-liquid corrections are 
taken into account. This makes thermal conductivity a 
very robust probe of the nodal quasiparticle spectrum in 
anisotropic superconductors. 

A detailed ARPES study on BSCCO [26] shows that 
Vf depends on angle 9 being maximum at 6 = 45° while 
kp only has a weak angle dependence (fcp=0.74 A -1 near 
9 = 45°). The vf(9) can be evaluated from a linear fit to 
the band dispersions (energy distribution curves) above a 
kink energy (about 50 meV below the Fermi level) . From 
the middle column of Fig. 4 in Ref. [26] , we estimate that 
hv F =l.4Q eVA for 9 = 45° and %v F =l.W eVA for 9 = 
38.5°. If we linearly extrapolate Vf with 9, we obtain 
Hvf = 1-09 eVA at an extended s-wave gap nodal angle 
of 9 n = 36.7° deduced above for the optimally doped 
top Cu02 layer of BSCCO. Similarly, we estimate hv F 
= 0.69 eVA for 9 = 25.8°, which is one of the nodal 
directions of the extended s-wave gap function deduced 
for slightly overdoped YBa2Cu307_ y (see below). It is 
interesting to note that the value of v F (1.05 xlO 5 m/s) at 
the nodal directions of YBa2Cu307_j, (YBCO) estimated 
from ARPES is very close to the value, (1.2±0.2) x 10 5 
m/s, estimated from the field dependence of the in-plane 
magnetic penetration depth at low temperatures [27]. 



With the values of hvp and hp, we can now calcu- 
late k /T for the d-wave and the extended s-wave gap 
functions using Eq. 3. The calculated values of k q /T are 
shown in the last column of Table 1. Since the top Cu02 
layer of this slightly overdoped BSCCO is optimally 
doped, as discussed above, we should compare the calcu- 
lated values of n /T with the measured one for an opti- 
mally doped BSCCO, which is not available. Fortunately, 
it is known that the value of k /T for slightly underdoped 
cuprates is slightly higher than that for slightly over- 
doped cuprates. This can be seen clearly from the YBCO 
system: k /T= 0.17±0.01 mW/K 2 cm for slightly under- 
doped YBa 2 Cu 3 6 .9o (Ref. [28]) and k q /T= 0.12±0.02 
mW/K 2 cm for overdoped YBa2Cu307.o [5]. By anal- 
ogy, we should take the experimental value of n /T to 
be about 0.2 mW/K 2 cm for optimally doped BSCCO, 
which is slightly larger than that for slightly overdoped 
BSCCO (0.15 mW/K 2 cm [5]) 

From Table 1, one can clearly see that the predicted 
value of k /T from the d-wave gap function is larger than 
the experimental value by a factor of 2 while this value 
from an extended s-wave gap model is within 10% of the 
experimental value. This indicates that the gap symme- 
try for the optimally doped BSCCO is not d-wave but 
extended s-wave. 

For optimally doped YBa2Cu30g.93, the measured 
k /T is about 0.17 mW/K 2 cm (Ref. [28]) and Am — 30 
meV. Taking the d-wave gap form of Eq. 2 with B = 0. 182 
and Am = 30 meV, we have n /T = 0.65 mW/K 2 cm, 
which is larger than the measured one by a factor of 4. 

On the other hand, an extended s-wave gap function 
A(0) = 24.5(cos46>+0.225) meV has been deduced from a 
single-particle tunneling spectrum for slightly overdoped 
YBCO [9]. This gap function has line nodes located at 
n = 25.8° and 64.2° in the first quadrant. With fiv F = 
0.69 eVA (see above), fc F = 0.74 A" 1 , and S = 47.7 meV 
(evaluated from the extended s-wave gap function), we 
calculate n /T = 0.12 mW/K 2 cm, in quantitative agree- 
ment with the measured one (0.14±0.02 mW/K 2 cm) [5]. 

The above quantitative data analyses unambiguously 
show that the gap symmetry for optimally doped 
cuprates is extended s-wave. How could this conclusion 
be compatible with all the phase-sensitive experiments? 
This issue can be resolved if we consider the fact that 



TABLE I. The calculated values of n a /T for the optimally doped BSCCO in terms of the d-wave and extended s-wave gap 
functions deduced from the best fits to the ARPES and FT-STS data for the optimally doped top layer of BSCCO (see text). 
The experimental value of k /T for optimally doped BSCCO is about 0.2 mW/K 2 cm (see text). Here d/n is the average 
separation between Cu02 planes stacked along the c axis, 9 n is the angle of the nodal direction, S = dA(9)/d8 is the slope at 
the node and N is the number of nodes per quadrant. 





d/n (A) 


9n 


S/N (meV) 


hv F (eVk) 


k F (A- 1 ) 


Vf/V2 


K /T (mW/K 2 cm) 


d-wave 


7.72 


45° 


21.38 


1.40 


0.74 


48.5 


0.378 


Extended s-wave 


7.72 


36.7° 


28.6 


1.09 


0.74 


28.2 


0.220 
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there are two types of charge carriers; one is intersite 
bipolarons of oxygen holes and another is Fermi-liquid 
type with a large Fermi surface [29]. The Fermi- liquid 
component is nearly absent for the hole doping p < 0.1, 
and increases monotonically with doping for p > 0.1 
(Ref. [29]). Further, it is shown that the Bose- Einstein 
condensate of bipolarons has d-wave symmetry [30] . Be- 
cause the interfaces of grain-boundary junctions consist 
of undcrdoped cuprates [31], the d-wave component of 
Bose-Einstein condensate of bipolarons is dominant at 
the surface, in agreement with phase-sensitive experi- 
ments based on grain-boundary Joscphson junctions [6]. 
It was also shown that the surface layer of a cuprate 
crystal is underdoped when it is contacted with a normal 
metal [32]. This can explain the observation of dom- 
inant d-wave component in the corner SQUID experi- 
ments [33,34]. The extended s-wave symmetry for the 
Fermi-liquid component can naturally account for phase- 
sensitive experiments based on out-of-plane Josephson 
tunneling [13]. 

In summary, we have analyzed data from angle 
resolved photoemission spectroscopy, Fourier trans- 
form scanning tunnelling spectroscopy, and low- 
temperature thermal conductivity for optimally doped 
Bi 2 Sr2CaCu208+ 2/ in order to discriminate between d- 
wave and extended s-wave pairing symmetry. The com- 
bined data are inconsistent with d-wave symmetry, but 
quantitatively consistent with extended s-wave symme- 
try with eight line nodes. 
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